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Abstract 

We discuss the history dependence and memory effects which are observed in the out- 
of-equihbrium conductivity of eiectron giasses. The experiments can be understood by 
assuming that the iocai density of states retains a memory of the sample history. We 
provide analytical arguments for the consistency of this assumption, and discuss the satu- 
ration of the memory effect with increasing gate voltage change. This picture is bolstered 
by numerical simulations at zero temperature, which moreover demonstrate the incom- 
pressibility of the Coulomb glass on short timescales. 

1 INTRODUCTION 

Classiness of localized electron systems was predicted long ago 00. Experimentally, such 
behaviour was first found in the slow capacitance relaxation of uncompensated semiconduc- 
tors 0, and more recently in a series of experiments on strongly disordered indium-oxide 
films 0, as well as in granular metallic films 5 where slow logarithmic relaxation, aging and 
memory were observed in the gate-controlled conductivity. 

The memory of equilibration conditions is manifested as an anomalous (symmetric) field 
effect. The conductivity of the sample always increases no matter whether the carrier density is 
increased or decreased. This out-of-equilibrium dip in the conductivity is remarkably universal, 
its shape being independent of disorder and magnetic field, while its width (in the applied 
gate voltage) shrinks with decreasing temperature and increases with carrier density A . The 
latter fact demonstrates the importance of electron-electron interactions, suggesting that these 
systems constitute genuine electron glasses. 

"X 

^ . 2 PHENOMENOLOGICAL THEORY OF THE MEMORY EFFECT 

_ _ i 

In Ref. 0, we proposed a phenomenological theory to account for these features. The theory is 
based on the standard model of Anderson insulators (in 2D) with strongly localized, classical 
electrons 

e 2 

H = 2_i n i~ v ) (nj - v) + y j nj(j)i, (2.1) 



the impurity sites being located on a square lattice with nearest neighbour distance a. Dis- 
order and frustration are introduced by random site energies (f>i , independently and uniformly 
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distributed over the range — W < <fii < W, and v is a homogeneously distributed background 
charge, adjustable by the gate. 

Three key ingredients are essential in order to understand the memory effect: 
i) Below a critical temperature (T < T g ) the electron glass falls out of equilibrium on 
experimental time scales due to the appearance of a large number of metastable states which 
are separated by high barriers. After a sudden change of external parameters the system is 
unable to relax rapidly to its new ground state, and may retain a memory of its former state. 
This scenario is suggested by numerical simulations [S], as well as by mean field theory 
which predicts a glass transition at T g ~ (e 2 / na) 2 /W for large disorder W 3> e 2 /na [TU] . 

ii) The repulsive interactions induce a Coulomb gap in the local density of states, which 
tends to the linear Efros-Shklovskii pseudogap at low temperatures [7j- The shape of the gap 
varies slightly from state to state and can serve as a fingerprint of the sample history. 

Hi) The system conducts via variable range hopping which is exponentially sensitive to 
perturbations of the local density of states. 

The theory of Ref. [fJJ makes the key assumption that the precise form of the Coulomb gap 
varies between different metastable states, and may thus reflect the pathway which a given 
state was obtained. This situation can occur when a small gate voltage is applied after a long 
equilibration. In the first place, the gate introduces additional background charges (changing 
v — > v + Av in (|2.1|0 . forcing the same amount of charge to enter the system in order to 
compensate on average for the gate field. The new carriers occupy the empty states just above 
the Fermi level, Ep, shifting the Fermi level up with respect to the bulk of the occupied states. 
The resulting density of states is asymmetric around the new Fermi level. We argue below 
that this feature survives the fast relaxations triggered by the presence of the new particles. 
Due to the asymmetry, thermal excitations of electrons into empty states will be enhanced, 
and the hopping conductivity will thus increase. A quantitative analysis shows that this effect 
dominates over the decrease due to the suppression of hole excitations [H]- 

3 METASTABILITY OF OUT-OF-EQUILIBRIUM STATES 

The introduction of new carriers renders other particles unstable and induces spontaneous 
hops that decrease the total energy. The question arises whether, or to what extent, such fast 
relaxation processes can wipe out the initial asymmetry of the density of states. In Ref. [Hj, 
we argued that as long as the density of new carriers, ny, is smaller than that of thermally 
excited electrons, ny, the gate voltage effect is perturbative. This condition is equivalent to the 
requirement that the shift of the Fermi level (with respect to the bulk of the states) be smaller 
than temperature. Below, we will provide more detailed arguments to justify the neglect of 
fast relaxations in this regime. On the other hand, at higher gate voltages new carriers are 
introduced on sites that were essentially always empty in the original state. Since in general 
their local environment is not favorable to the addition of a particle, the new electrons likely 
trigger relaxation processes and destabilize the original state. Such reconfigurations will lead 
to the saturation of the anomalous field effect at voltages where ny > nr- 

Let us analyse the stability criterion ny < np in more detail in the case of low tempera- 
tures where the equilibrium density of states, po(£ = E — Ep), displays a linear Coulomb gap. 
More precisely, we assume the scaling form po(£) = T/{e 2 / n) 2 f{\£\/T) with /(0) = 0(1), and 
f(x) ~ x for 1 <C x <C (e 2 / ' na) 2 jWT '. If the density of injected particles, ny, is larger than 
rip ~ (kT /e 2 ) 2 , their mutual interaction exceeds the temperature, so that they spontaneously 
rearrange, reshuffling the distribution of levels. On the other hand, for ny < np, the interac- 
tions between new particles can be neglected. The particles will a priori occupy the first empty 
states above Ep, shifting the Fermi level by AEp ss Av/po(0) ~ (e 2 / n) 2 Ap /T , so that the 
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new density initially takes the form pi(£) — po(£ + AEp). However, the Coulomb repulsion 
from the new electrons shifts the local energies of other sites, inducing spontaneous hops that 
affect the distribution of levels. 

Let us estimate the density n p (E*) of recombining electron- hole pairs that affect the number 
of states in the interval < |£| < E*. We count the pairs (ij) whose recombination energy, 
AEij = e 2 / Hir^ — — is raised above temperature, while the pair was stable (AE^ < T) 
before the introduction of new particles. We approximate n p (E*) by the density of states for 
pairs with AEij of order T, multiplied by the typical shift of A-Ey for a pair of sites with 
distance ry. The latter is of order nj Ave 2 j k for the dominant short pairs with nj < Av^ 1 ! 2 . 
Assuming the single-site energies to be independent, we find 

n p (E*) « J d 2 r J* d£ipi(£i) J° ^ d£ 2Pl (&)rAv^5 - & - \&\-T\ . (3.1) 

This should be compared to the change nd(E*) in the density of particle/hole states in the 

same energy interval |£| < E*, caused by the shift of E F : n d (E*) s» (pi(£) — pa(£))dt; ~ 
E* AEp I (e 2 j k) 2 . Performing the integration over r in (|3.1I) . and dividing by rid(E*), we obtain 

n p (E*) n f E * f° 1 e 4 Av T 

Me^~Wae^J ^M a) j E , (fl + 16,1 + r)» fx + 16| + r ~ £f (3 ' 2) 

which is only a small fraction in the range of energies (E* 3> T) that are probed in hopping 
conductivity. Note, however, that within the range |£| < T the density of states can be 
appreciably modified by fast relaxations. 



4 NUMERICAL SIMULATIONS AT T = 



In order to provide evidence that history dependence can be seen in the density of states, we 
carried out simulations at T — for the model (|2.1|l with moderate disorder strength W = 0.6 
in systems of linear size L = 32, 48, 64, 80. Starting from random initial conditions, we allowed 
all hops that decrease the energy, until we reached a state that is stable to any single hop. We 
then filled in a fraction Av <C v « 0.54 of new particles, putting them onto the available sites 
with least energy cost. Finally, we simulated the fast relaxation processes by steepest descent 
via single hops to the closest metastable state. For both metastable states we determine the 
Fermi level Ep (the energy separating occupied and empty levels) and calculate the distribution 
of local energies, Ei — dH/drii. For the states obtained before particle injection, Ep is found 
to increase linearly with filling fraction, dEp jdv « 2.45, as one expects for a system with finite 
("field cooled") compressibility. However, in the states obtained after injection and relaxation, 
the average Fermi level is systematically higher than that of unbiased states with the same 
filling fraction. This is also manifested by the fact that the bulk of the density of states, pi{£), 
is shifted down by AEp, see Fig. ^ The extra shift AEp increases with filling fraction as 
Ai/ 1 / 2 , as one expects from the presence of a linear Coulomb gap. This result reflects the 
incompressibility of the Coulomb glass 0, k = (dAEp / dAv)^ 1 ~ {Av) 1 / 2 — * (for small 
Av). A similar phenomenon occurs in the pinned Wigner crystal At larger Av, the shift 
saturates to AEp°, as the system becomes unstable and reorganizes significantly. 

At finite, but small temperature the numerical results for T — should remain robust, since 
the density of states can only relax further by activated processes. The fact that at T = the 
shift AEp saturates only at a finite fraction Av suggests that the saturation condition ny ps ut 
for the anomalous field effect only applies for T > AEp°. However, at least in the presence of 
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Fig. 1. Left: The density of states before and after particle injection, plotted as a function of 
£ = E — Ef- The data are averaged over 200 samples with 80x80 sites, for disorder with width 
W — 0.6. Right: The shift AEf as a function of the relative increase of carrier density. AEf is 
determined from the shift between the bulk of the density of states. It agrees within error bars with 
the difference between the average Fermi level in both kinds of states with the same filling fraction. 
The full line is the best fit to the form AEf = a^^Au/u)- 1 + c" 1 )" 1 . 

strong disorder, W 3> e 2 /na, one expects a rather large temperature interval AE^° < T < T g 
in which the width of the anomalous field effect is correctly described by ny ~ tit ■ 

In conclusion, we have provided quantitative arguments and numerical evidence at T = 
that the local density of states can exhibit memory of the sample history. These independent 
approaches provide further support for the memory scenario proposed in Ref. 
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